Wave packet revivals in a graphene quantum dot in a perpendicular magnetic field 
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We study the time-evolution of localized wavepackets in graphene quantum dots in a perpendicular 
magnetic field, focusing on the quasiclassical and revival periodicities, for different values of the 
magnetic field intensities in a theoretical framework. We have considered contributions of the two 
inequivalent points in the Brillouin zone. The revival time has been found as an observable that 
shows the break valley degeneracy. 
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I. INTRODUCTION 

Graphene, a single sheet of graphite, has attracted 
growing interest due to its remarkable and starling prop- 
erties and its potential applications in nanoelectronics 
P-Q. Many of these properties are attributed to the 
existence of a peculiar band structure constituted by 
quasifree electrons being described by a massless Dirac 
equation Q. In fact, perfect tunneling of electron wave 
packets through potential barriers (the Klein paradox) 
Q or the unconventional quantum Hall effect, have been 
experimentally confirmed in graphene Q. It has been 
also argued that graphene in a magnetic field is a promis- 
ing system for the experimental observation of the highly 
oscillatory motion of the electric current known as zitter- 
bewegung [§-[l2j ■ In addition, revivals of electric currents 
in graphene in the presence of a magnetic field [l3j have 
been theoretically predicted and an exhaustive study of 
quantum wave-packet revivals, fractional revivals, classi- 
cal periodicity in graphene was reported (l3l . fl4j . 

In recent times, graphene quantum dots have been 
widely studied, both theoretically and experimentally. 
It has been reported that there are different theoreti- 
cal approaches which allow to electrostatically define a 
graphene quantum dot, or to create it by cutting flakes 
of graphene 15 22]. Also, there has been an analysis 
of the level statistics of the graphene quantum dots ob- 
served in experiments (for quantum dots smaller than 
100 nm) which has concluded that it is well described by 
the theory of chaotic Dirac neutrino billiards [23[ . On the 
other hand, new experiments have been designed for the 
creation of graphene quantum dots (23 - [26| by etching or 
scratching graphene islands including the study of energy 
levels in magnetic fields (2(| [27j ■ Additionally, graphene 
quantum dots have been point out to be very attractive 
as spin qubits in quantum information processing [28]. 

The temporal evolution of wavepackets, relativistic and 
nonrelativistic, displays interesting revival phenomena 
due to quantum interference. Several types of periodic- 
ity may emerge depending on the energy eigenvalue spec- 



trum, in this regard revivals and fractional revivals have 
raised great interest during recent years [2{|. Propagat- 
ing wave packets initially evolve quasiclassically and os- 
cillate with a period T c i , and later the wavepackets spread 
and collapse. For longer times, the wavepackets evolve 
in the so-called revival time T rev to a state that closely 
reproduces its initial waveform. Fractional revivals ap- 
pear as a temporal self-splitting of the initial wavepacket 
into a collection of a number of scaled copies at times 
t = pT rev /q, with p and q mutually prime. Assuming that 
the initial wavepacket is a superposition of the eigenstates 
u n sharply peaked around some no, revival times can be 
obtained from the Taylor series of the energy spectrum 
E n around E na as T rev = AttH/\E'^ |. Revivals have been 
investigated theoretically and experimentally in nonlin- 
ear quantum systems [29l - l34| . Recently, revivals of elec- 
tric current in monolayer graphene in the presence of an 
external magnetic field have been also described [l3|, [TiJ . 

This article is organized as follows. In section II the 
phenomenon of wavepacket revivals in a graphene quan- 
tum dot in a magnetic field using a model Hamiltonian is 
studied. We discuss on the relevance of quantum revivals 
for the characterization of the breaking of the valley de- 
generacy in graphene quantum dots. 

II. WAVE PACKET REVIVALS FOR 
GRAPHENE QUANTUM DOTS IN A 
MAGNETIC FIELD 

To study wave packet revivals in a graphene quantum 
dot in a perpendicular magnetic field B we consider the 
model Hamiltonian for electrons in the valley-isotropic 
form which is given by (2l| 

H T =v F (p + eA)-cr + TV(r)a z . (1) 

We use the symmetric gauge for the vector potential, 
A = B/2(— y, x, 0) = B/2(— r sin </>, r cos <j>, 0), where <j> 
the polar angle, vf = 10 6 m/s the Fermi velocity, and 
t = ±1 differentiates the two valleys K\ and K%. <x are 
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FIG. 1. Positive energy spectrum for a graphene Quantum 
Dot, with R = 70 nm, versus the quantum number m and 
different values of the magnetic field B (Ki point). 



B{T) 


T cl (ps) 


T rev (ns) 


r d (Ps) % ev (ns) 


0.1 


0.14 


115.8 


0.14 


120.7 


0.5 


0.14 


6.8 


0.14 


7.6 


1 


0.14 


2.3 


0.14 


2.6 


4 


0.14 


0.2 


0.14 


0.3 



TABLE I. Classical period T c ; and revival time T rev in the 
two valleys (Ki and K2) in the band structure of graphene 
quantum dot for different values of the magnetic field B and 
for an initial wave packet with no = 20, m = 0, a = 1.1. 



Pauli's spin matrices in the basis of the two sublattices of 
A and B atoms. We assume the confinement potential in 
21], which is a mass- related potential energy V(r) cou- 
pled to the Hamiltonian via the a z Pauli matrix, where 
V(r) — for r < R and V(r) = 00 for r > R, that is, 
tends to infinity in the edge of the dot. 
Upon introducing the magnetic length Ib 



and using the fact that H T commutes with the total an- 
gular momentum operator J z = L z + |er z , [H T , J z ] = 0, 
the solution of the Dirac equation H T ip(r, <fi) — Eip(r, </>) 
(where ip(r,4>) = h/>i(r, <j>), fair, <ft)1 is a two-component 
spinor) is given by 21.]: 
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is a normalization constant. The characteristic equation 
for the allowed eigenenergies E of the Quantum Dot u sing 
the boundary condition for a circular confinement (2ll . |35| 
1P2/1P1 — Texp[i<fi] has been obtained in j2l| and it can 
be written as 
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FIG. 2. Time dependence of |yl(t)| 2 for an initial spherically 
symmetric Gaussian wavepacket with no = 20, a = 1.1, B = 
0.2 T and R = 70 nm for (a) the valley Kx (T rev = 28 ns) 
and (b) the valley K2 {Trev — 31 ns). 



The spectrum of electrons confined in a circular 
graphene quantum dot in a perpendicular magnetic field 
is given by E(n, m, r) verifying @ and ^ and the en- 



(to + l),m,r 2 /2lj 



ergy eigenfunctions are ip( n ,m)( r i 4 1 )- Schnez et al. [21 1 
have analyzed the limits for B = and R/Ib 00 cases. 
For B = the eigenenergies are given by the charac- 
teristic equation rJ m (kR) = J m+ \{kR). On the other 
hand, they showed that the energy leve ls converge to the 
landau levels E m — hvFk m = ±VF\/2ehB(m + 1) for 
R/Ib 00 ■ 

The initial wave packets will be constructed as a super- 
position of eigenstates assuming a Gaussian population 
of energy levels 



x Zll^L ( ^ - (to + 1), m + 1, r 2 /2l% 



(4) 



kl B V 2 

(2) that is, we choose c„ i7n Gaussianly distributed and cen- 
where L(a, b, z) is the generalized Laguerre function and c tered around some no, Too with width a (See for instance 



3 



Refs. [33, |36j, |37[ ) • We will study wave packet dynamic of 
this initial wave packet for different values of the inten- 
sity of the magnetic field B. The particular assumption 
we have made to construct the initial wavepacket has 
an influence in the results (see i. e. [32J). Moreover, the 
construction of a realistic initial wavepacket (i.e. by laser 
excitation) would require further considerations. 

We start studying the special case of a radially sym- 
metric wave packet, that is, a superposition of purely 
to = states. For an initial wave packet sharply peaked 
around a central no, the different time scales can then be 
identified from the coefficients of the Taylor expansion of 
the energy spectrum E n around E no : 
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In Table I (middle column) the times scales T c i and T rev 
have been numerically fitted from energies of Fig. f. The 
quantum number dependence of the energy eigenvalues 
is necessary to determine the classical period and the 
revival time in a system. In Fig. Q]we plot the n and m 
dependence of the energies E n ^ m for B = 0.1, 0.5, 1 and 
4 T and for R — 70 nm in the valley K\ . 

To analyze the evolution of the wave packet we shall 
use the autocorrelation function 



A{t) 



r dr d<p ip*(r, 4>, Q)ij)(r, <j), t). 



(6) 



An alternative approach based on information entropies 
has been recently proposed [38|. The occurrence of re- 
vivals correspond to the approximate return of |^4(£)| 2 
to its initial value of unity. The time evolution of the 
autocorrelation function for an initial wave packet with 
n = 20, to = a = 1.1, R = 70 nm and B = 0.2 T is 
showed in Fig. [5] (top panel). We can see that at early 
times the Gaussian wave packet evolves quasiclassically 
in the inset of the figures. For longer time scales the wave 
packet initially spreads and delocalizes. A revival of the 
wavepacket can be observed for T rev /2 (mirror revival 

Let us turn to examine quantum revivals in the two 
inequivalent points in the Brillouin zone K\ and Ki 
which corresponds with Hamiltonian (TT]) with r = 1 and 
t = — 1 , respectively. We have calculated the classical 
and revival times for different values of the magnetic field 
B which are summarized in table HI The quasiclassical 
evolution is analogous in both valleys as we can see in 
the insets of Fig. [5J We have found both theoretically 
and by simple inspection of the behavior of the auto- 
correlation function in Fig. [5J that T rev differs in both 
zones clearly for a magnetic field B = 0.2 T. In fact, 
it is known that the valley degeneracy is broken by the 
perpendicular magnetic field [20j,|28[. So, we have found 
the revival time as an observable which shows explicitly 
the breaking of valley degeneracy in a graphene quan- 
tum dots. We would like to point out that the breaking 
of the valley degeneracy is a prerequisite for spin-based 
quantum computing in graphene quantum dots . As the 




FIG. 3. Time dependence of |j4(£)| 2 for initial wavepackets 
with B = 4 T, R = 70 nm, no = 1, mo = —15 and a n — a m = 
0.7 in the valley Ki.(a) First classical periods of motion with 
T c ( "' = 0.08 ps, = 0.65 ps (b) Long time dependence 

with T, ( ™2 = 2.09 ps and TrTv = 14.50 ps. 



magnetic field is growing the revivals times are equalizing 
as we expect for Landau levels regime 

Finally, we consider an initial superposition state with 
coefficients taken from a Gaussian distribution centered 
around to = — 15 and n = 1 for B = AT. In this case 
the energy levels are characterized by two quantum num- 
bers E(n,m). Systems with two quantum numbers offer 
richer possibilities for wave packet revivals (for a the- 
oretical analysis see [32[ ). The classical periods yield 
straightforwardly T£ = 0.078 ps and T c 7 = 0.65 ps (from 
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Now, the long-time revival structure de- 
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and the revival struc- 



AttH/ 

T^ m) = 2nh/ 

ture is constructed on the interplay between these times. 
This is shown in Fig. [3] for B — AT. In the top panel 
it is apparent that at short times the quasiclassical os- 
cillatory behavior with period T c j = 0.078 ps. At 
medium times the quasiclassical behavior corresponding 

to T^™ 1 ' = 0.65 ps is superimposed on the first oscil- 
lations. In the bottom panel revivals can be seen at 
t w 8.5ps as an interplay between the revival times, as 
we would expect. Snapshots of the numerical simulation 
of the probability density are given in Fig. [4] at several 
times. It is apparent that the regeneration of the density 
function is given at time t = 8.65 ps. 

We want to remark that the typical scattering times 
should be taken into account for a more realistic descrip- 
tion of this problem. Moreover, it would be interesting 
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FIG. 4. Snapshots of the probability density \ip{x, y)\ 2 for an 
initial gaussian wavepacket with no = 1, mo = —15, a n = 
a m — 0.7, B — 4 T and R — 70 nm. It is noticeable the 
existence of fractional revivals after the derealization of the 
initial wave packet (middle panels) and the revival at t = 
8.65 ps. 



to understand the influence of shape deformations, im- 
purities, or interactions on revivals, so a more detailed 
study of this topic could be made in a analogous way 
that has been done before (see i.e. [l3|, [TJ]) in graphene 
in a perpendicular magnetic field. 



III. CONCLUSIONS 

In this work we have presented time-dependent effects 
in the propagation of wavepackets on graphene quantum 
dots for different values of the magnetic field. When the 
wavepacket is sufficiently localized around some central 
quantum numbers (no, Too), it evolves quasiclassically. 
At later times, it spreads leading to the collapse of the 
classical oscillations, but at times that are multiple or ra- 
tional fractions of T rev , the wavepacket regains its initial 
form. We have also shown revival time as an observable 
sign of breaking the valley degeneracy, that it is known 
as a prerequisite for spin-based quantum computing. It 
can be underline that we have found revivals of the elec- 
tron wavepacket in a Dirac system and we suggest that 
this realistic system could be used to simulate relativistic 
quantum revivals [3, Eil . 
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